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Abstract 

We  give  results  of  recent  work  on  a  newly  developed  frequency  stability  characterization, 
called  Total  variance,  whose  main  advantages  are  improved  confidence  at  and  near  the  longest 
averaging  time  of  half  the  data  duration,  and  lower  sensitivity  to  drift  removal.  Properties  given 
here,  for  the  standard  FM  noise  types,  include  mean,  degrees  of  freedom,  frequency  response, 
and  empirical  distribution  function. 


1  Introduction  and  Conclusions 

This  paper  is  about  characterizing  common,  difficult  to  characterize  frequency  noise  modulations 
found  at  long-term  averaging  time  r  in  the  output  signal  of  many  laboratory  frequency  standards. 
It  assumes  a  familiarity  with  the  Allan  variance  and  its  characterizations  of  white,  fficker,  and 
random  walk  FM  noise  models  (WHFM,  FLFM,  and  RWFM)  [1]. 

A  shortfall  of  the  currently  recommended  [2]  Allan  variance  for  r,  denoted  here  by  Avar  (r), 
is  that  the  usual  estimators  of  it  are  highly  variable  at  large  r  [1,  3,  4]  and  are  sensitive  to  the 
method  of  drift  removal  [5].  We  report  statistical  properties  of  a  new  kind  of  frequency  variance, 
inherently  dependent  on  measurement  duration  T  as  well  as  r,  called  Total  variance  and  denoted 
hereafter  by  Totvar  (r,  T)  (pronounced  tot  '-var).  We  quantify  the  improvement  in  the  uncertainty 
on  frequency-stabihty  estimation  from  the  use  of  Totvar  (r,  T)  rather  than  Avar  (r)  in  the  presence 
of  FM  noises.  The  square  root  of  this  variance,  cahed  Total  deviation  and  denoted  by  Totdev  (r,  T) 
or  the  recommended  [2]  notation  total  ("t),  can  be  interpreted  like  the  Allan  deviation  <Tj,(t), 
but  with  improved  confidence  at  long-term  r,  as  pointed  out  in  earlier  papers  [6,  7]  and  quantified 
in  this  paper  for  FM  noises.  We  do  not  address  the  important  property  that  Totvar  (r,  T)  appears 
to  have  considerably  less  sensitivity  to  the  method  of  drift  removal  than  Avar(r)  [8]. 

The  main  advantages  of  Totvar  over  Avar  are  improved  confidence  at  and  near  the  longest 
averaging  time  of  r  =  r/2,  and  lower  sensitivity  to  drift  removal.  By  theory  and  simulation  we 
have  computed  its  mean,  variance,  and  empirical  probability  distribution  in  the  presence  of  the 
three  FM  noise  types.  Variance  results  are  given  in  terms  of  equivalent  degrees  of  freedom  [9].  In 
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the  presence  of  white  FM  noise  modulation,  Total  variance  is  an  unbiased  estimate  of  the  Allan 
variance  for  all  r,  yet  has  three  degrees  of  freedcon  instead  of  one  at  r  =  T/2.  For  all  three  noise 
types,  the  mean  and  edf  of  Totvar  are  given  by  simple  exact  or  empirical  formulas.  A  comparison  of 
the  empirical  distribution  of  Totvar(T/2,T)  with  the  appropriate  chi-squared  distribution  indicates 
that  confidence  intervals  based  on  chi-squared  levels  are  conservative. 

In  the  estabhshed  tradition  of  time  and  frequency  statistics,  frequency  stability  is  described 
in  terms  of  finite-dilference  variances  that  are  ensemble  or  infinite-time  averages  of  stationary, 
ergodic  increments  of  phase  [12,  13,  14].  In  particular,  the  theoretical  Allan  variance  is  a  number 
that  depends  only  on  r,  while  its  conventional  estimators  are  random  variables  that  depend  on 
both  T  and  T,  the  largest  possible  r  being  T/2.  Total  variance  is  a  random  variable  that,  along 
with  all  its  properties,  inherently  depends  on  both  r  and  T,  Moreover,  Total  variance  can  report 
values  beyond  the  usual  Allan  variance  last-r  value  of  r  =  T/2;  its  values  at  r  >  T/2  might  be 
used  to  augment  the  normal  last-r  value  of  frequency  stability  reported  at  r  =  T/2, 

We  compute  the  frequency  response  of  Total  variance  as  a  function  of  r  by  averaging  the 
squares  of  the  Fourier  transforms  of  Total  variance  sampUng  functions,  and  find  that  it  resembles 
the  frequency  response  of  Allan  variance.  The  results  of  these  investigations  indicate  that  Total 
variance,  while  it  has  an  interpretation  like  that  of  the  Allan  variance,  also  has  lower  variability 
and  less  sensitivity  to  drift  removal. 

2  Equations  for  Totvar  (r,  T) 

The  purpose  of  this  section  is  to  give  a  precise  definition  of  Totvar  (r,  T)  for  an  iV^-point  time 
deviation  record  with  sample  period  tj).  In  the  following  description,  the  indices  m,  n,  and  Nx  are 
related  to  time  by  r  =  mro,  t  =  to  +  nro,  and  T  =  N^tq,  where  to  is  the  time  origin  and  without 
loss  may  be  made  equal  to  0. 

We  start  with  time-deviation  data  n  =  I  to  Nx,  with  normalized  frequency  deviations 
Vn  =  (iCn4-i  —  ^n)  /tq,  Ti  ~  1  to  Ny  =  Nx  —  1.  Extend  the  sequence  {j/n}  to  a  new,  longer  virtual 

sequence  {j/*}  by  reflection  as  follows:  forn  =  1  to  Ny  let  Vn  =  yn,  for  j  =  1  to  IVy  —  1  let 

Vi-j  =  yj,  y*Ny+j  =  VNy+i-j-  (1) 

An  equivalent  operation  can  be  performed  on  the  ori^al  time-deviation  sequence  {xn}  to  produce 
an  extended  virtual  sequence  {x*  }  as  follows:  forn  =  1  to  Nx  let  x*  =  Xn,  for  j  =  1  to  Nx  —  2  let 

Xi_j  =  2xi  Xi^j,  ^Nx+j  ~  ^^Nx  ^Nx~j‘  (2) 

This  operation,  depicted  in  Figure  1,  is  called  extension  by  reflection  about  both  endpoints.  The 
result  of  this  extension  is  a  virtual  data  sequence  x^,  n  —  3  —  JVi  to  2Nx  —  2,  having  length  3Nx — 4 
and  satisfying  y*  =  (x*^i  -  x* )  /tq,  n  =  3-Nx  to  2Nx  -  3, 

We  now  define 

j  Nx-l 

Totvar  (m,  Nx,  tq)  =  ^  “  2x*  -f  x*^J^ ,  (3) 

2  {niTb)  [Nx  ~  2)  „=2 

for  1  <  771  <  A*  -  1,  that  is,  r  can  go  to  (AT*  -  1)to  instead  of  the  usual  Umit  of  [(Ax  -  1)  /2J  tq. 
The  previous  notation  Totvar  (r,  T)  is  to  be  regarded  as  equivalent  to  (3)  with  the  dependence  on 
7t)  suppressed.  Totvar  can  also  be  represented  in  terms  of  extended  fractional  frequency  fluctuation 
averages  as 

Totvar  (m,  Ny,  tq)  =  Z)  H  “  Vn-m  . 
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(4) 


where  Tn  (jn)  =  -  x^)  /  (mro). 

It  can  be  verified  that  Totvar,  like  Avar  and  its  estimators,  is  invariant  to  an  overall  shift 
in  phase  and  frequency;  that  is,  if  a  first*degree  polynomial  co  +  cin  is  added  to  the  original  data 
set  Xt»,  then  Totvar  does  not  change. 

We  expect  the  Total  variance  to  be  applied  mostly  to  long  phase  records  in  which  the 
FM  noises  dominate  the  PM  noises  (white  and  flicker  PM).  Then,  it  is  convenient  for  theoretical 
purposes  to  approximate  Totvar  (m,A/a;,ro)  by  a  continuous- time  analog  called  Totvarc(r,r),  in 
which  the  sum  in  (3)  is  replaced  by  an  integral  and  dependence  on  to  is  eliminated.  The  time 
deviation  is  now  a  continuous-time  process  x{t),  given  for  0  <  t  <  T.  Extend  x(t)  to  x*{t) 
as  follows:  for  0  <  t  <  T  let  x*  (t)  =  x(t);  for  0  <  it  <  T  let  x*  (—it)  =  2x(0)  —  x(u)  and 
x*{T  +  u)  =  2x  (T)  “  X  (r  -  u).  Then  x*  (t)  is  given  for  -T  <  t  <  2T.  The  continuous  analog  of 
Totvar  is  defined  by 

Totvarc  (r,  T)  -  ^  [x*  (t  ~  r)  -  2x*  (t)  +  x*  (t  -h  r)]^  dt,  (5) 


for  0  <  r  <  r. 

The  expressions  above  are  quite  different  from  their  equivalent  Allan  variance  expressions. 
Mainly,  Totvar  (r,  T)  reports  a  value  for  an  interval  r  within  data  length  T  based  on  more  samples 
of  the  second-difference  of  phase  (or  first-difference  of  average  frequency)  using  a  rearranged  and 
extended  series  of  the  original  data  series  {xn}*  It  does  this  by  a  multiple  sample  on  the  phase 
using  a  larger,  virtual  set  of  data  which  is  an  odd,  or  reflected,  extension  at  the  beginning  and  end 
(left  and  right)  of  the  original  real  set.  Figure  1  illustrates  the  extension  and,  hence,  the  resulting 
circular  or  repeating  representation. 

Analyzing  a  larger  virtual  data  set  built  from  the  original  data  set  has  been  a  tool  in 
frequency-domain  signal  processing  for  many  years.  An  assumption  of  periodicity  replaces  the 
recurrent  behavior  (in  a  time-series  sense),  a  consequence  of  the  ergodic  principle  [15].  In  particular, 
a  range  of  frequency  values  (Fourier  components)  from  0  <  f  <  can  be  extended  by  a  mirror 
reflection  through  /  =  0  so  that  “negative  frequencies”  are  added  to  an  original  data  set,  resulting 
in  -  fk  <  f  <fh-  In  the  context  of  time-series  analysis,  rather  than  doing  extensions  of  the  original 
vector  {xt,}  and  applying  the  straight  second-difference,  we  alternatively  can  resample  within  the 
original  vector;  see  Section  4  for  an  algorithm  that  requires  no  extension  of  the  data  array. 


3  Totvarc  (r,  T)  as  an  Estimator  of  Avar  (r) 

3.1  Mean  and  Variance 


For  computing  theoretical  moments  it  is  convenient  to  use  the  continuous-time  random  variable 
Totvarc  (TjT)  as  a  surrogate  for  the  discrete- time  random  variable  Totvar  (t,T),  for  the  same 
reason  that  the  calculus  of  integrals  and  derivatives  is  less  intricate  than  the  calculus  of  sums  and 
differences.  The  mean  and  variance  of  Totvarc  (r,  T)  in  the  presence  of  the  three  standard  power-law 
FM  noises  were  computed  by  the  generalized  autocovariance  method  [4]  under  the  assumption  of 
Gaussian,  mean-zero  second  differences  of  phase;  no  frequency  drift  or  drift  removal  is  allowed.  The 
mean  E  [Totvarc  (r,  T)]  is  compared  to  Avar(r);  the  variance  is  most  conveniently  communicated 
through  the  equivalent  degrees  of  freedom  (edf),  defined  for  a  random  variable  V  by 


edf(l^)  = 


2{EVf 
var(V')  ’ 


(6) 
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The  results,  some  of  which  are  exact  and  some  of  which  are  empirical  fits  to  numerically  computed 
results,  can  be  expressed  as 


E  [Totvarc  (r,  T)]  r 

Avar(r)  “  “t’ 

T 

0<T<-, 

(7) 

T 

edf  [Totvarc  (r,  T)|  =  b - c, 

T 

T 

0<r<-, 

(8) 

where  a,  6,  and  c  are  given  in  Table  I.  These  results  were  checked  by  simulations  of  Totvax  (m,  Nx,  tq), 
with  Nx  ~  101. 


Tbble  I.  Values  of  o,  6  and  c  for  FM  noises 


Noise 

a 

b 

c 

WHFM 

0 

3/2 

0 

FLFM 

1/(3  M2)) 

24  (ln(2)/7r)^ 

0.222 

RWFM 

3/4 

140/151 

0.358 

The  simple,  exact  form  (7)  for  the  mean  of  Totvar  can  be  interpreted  as  a  scaling  property 
of  power-law  noise.  It  turns  out  this  way  because  the  shapes  of  the  sample  functions  for  Totvar  and 
Avar  (see  Section  4)  depend  only  on  t.  For  Tf2  <t  <T  the  sampling  function  shapes  depend  also 
on  T;  yet,  it  is  noteworthy  that  (7)  persists  all  the  way  to  T,  but  only  for  white  FM  and  random 
walk  FM.  The  edf  results  are  empirical,  with  an  observed  error  below  1.2%  of  numerically  computed 
values;  for  white  FM,  though,  the  edf  result  appears  to  be  exact,  although  this  is  imverified. 

For  white  FM  noise  processes.  Total  Variance  is  an  unbiased  estimate  of  the  traditional  Allan 
Variance  (square  of  frequency  change  vs.  r)  for  all  averaging  times  (r’s).  Its  primary  advantage,as 
surmised  from  (8)  and  Table  I, is  a  considerable  improvement  in  the  confidence  of  that  estimate  at 
longer  averaging  times.  For  example,  10  000  seconds  of  frequency  measurements  means  that  r  can’t 
go  beyond  5000  seconds  (T/2)  to  get  one  single  estimate  of  frequency-change  over  the  data  duration 
(last  half  minus  first  half)  using  Avar.  The  edf  result  (8)  for  white  FM  yields  the  equivalent  of 
three  independent  estimates  (edf  =  3)  by  using  Totvar  and  its  multi-sampling  function,  which  is 
discussed  in  the  next  section.  The  improvement  of  edf  in  the  presence  of  FLFM  and  RWFM  (2.097 
and  1.514,  respectively)  is  not  as  dramatic,  but  substantial  nevertheless.  For  r  —  T/4,  T/8,  etc.,  the 
confidence  measures  of  Allan  and  Total  confidence  approach  each  other,  until  they  are  essentially 
the  same  at  T/16.  For  small  values  of  t/T,  the  time-deviation  record  is  extended  only  a  short 
distance  at  both  ends,  and  hence  Totvar  (t,  T)  differs  little  from  the  fully  overlapped  estimator  of 
Avar(T)  [9]. 

3.2  Distribution  Functions 

In  the  tradition  of  time  and  frequency  statistics,  it  is  customary  to  derive  confidence  intervals  for 
frequency  stability  on  the  basis  of  the  assumption  that  the  probability  distribution  of  a  frequency 
stability  estimator  V,  when  scaled  appropriately,  follows  the  chi-squared  distribution  with  the  same 
edf  as  V;  see  [9].  The  chi-squared  assumption  has  been  investigated  in  a  limited  way  for  conventional 
estimators  of  Allan  variance  [10,  11];  in  view  of  the  greater  complexity  of  Total  variance,  some 
investigation  of  its  distribution  is  appropriate. 

Let  V  denote  the  Totvar  estimator  of  Allan  variance  tr^,  for  some  r,  and  let  r  =  F  (V)  /a^, 
=  edf  (V),  which  are  presumed  known  from  the  previous  results.  Then  the  random  variable 


has  the  same  mean  and  edf  as  does  a  Xut  namely,  u.  If  X  had  the  distribution,  then  one  could 
derive  confidence  intervals  for  based  on  one  observation  of  V. 

We  carried  out  a  brief  investigation  of  the  chi-squared  assumption  for  the  most  important 
case  of  r  —  T/2.  One  thousand  independent  trials  of  Totvar  (50, 101)  were  simulated  for  the  three 
FM  noise  types,  and  the  empirical  distribution  functions  of  X  were  plotted  along  with  the  theoretical 
chi-squared  distribution  functions.  The  chi-squared  distributions  were  good  fits  to  the  empirical 
distributions  except  at  the  lower  tails,  expanded  views  of  which  are  shown  in  Figure  2.  The  empirical 
tail  always  lies  below  the  chi-squared  tail.  Thus,  for  a  probability  p  <  0.2,  if  (p)  is  the  chi-squared 
level  for  p,  and  x*  (p)  the  level  of  the  scaled  Totvar  X,  we  have  x^  (p)  <  x*  (p).  Because  the  upper 
end  of  a  true  confidence  interval  for  cr^  is  proportional  to  1/xl  (p)  for  an  appropriate  value  of  p, 
using  X  (p)  for  this  purpose  gives  a  conservative  confidence  interval.  For  example,  suppose  we  want 
a  90%  confidence  interval  for  <r^  at  r  =  T/2  based  on  F  =  Totvar  (T/2,  T)  and  a  white-FM  noise 
assumption.  Then  r  =  1,  =  3;  the  5%  and  95%  xi  levels  are  0.352  and  7.81;  and  a  90%  confidence 

interval  for  based  on  chi-squared  is  [3F/7.81,3V/.352]  =  [0.384V,  8.52  V].  (Take  square  roots  for 
(7.)  The  more  realistic  value  of  0.60  for  the  5%  Totvar  level  from  Figure  2(a)  reduces  the  upper  end 
of  the  confidence  interval  to  5V. 

4  Sampling  Functions  and  Frequency  Responses 

Extending  by  reflection  an  original  {x„}  vector  at  both  ends  and  then  applying  a  second  difference 
can  be  equivalently  represented  as  four  different  types  of  differencing  on  {x„  :  n  — 
directly.  The  summand  of  Totvar  (m,iVa:,To)  in  the  following  equation  takes  on  four  forms  that 
depend  on  the  relationship  of  n  to  a  given  m  and  Nxi 


Totvar  (m,  Nx,ro) 


1 

2(mTonNx-2) 


Nx-l 

E^n, 


n=2 


where 


Dn  —  2Xji 

”  2xi  X2— n+m  ^  2x,i  ■{“  Xn+m 

Dn  ~  2^71—77,  ~  2x71  V  2x^j  X2fiFx—n—m 

Dn  —  2X\  X2— 7i-(-m  2X7,  2xfj^  ^2Nx—n—m 


m  —  n>l,  m  +  n  <  Nx, 
m  —  n<l,  m  +  n  <  Nx, 
m  — n>l,  m  +  n>  Nx, 
m  —  n<l,  m  +  n>Nx. 


(9) 


We  can  derive  ftom  these  expressions  the  frequency  sampling  functions  associated  with 
Totvar (r,T),  that  is,  how  its  terms  act  on  Pn  =  (x,,  —  Xn-i) /tq,  and  contrast  them  with  the 
simpler  sampling  function  associated  with  Avar  (r)  (see  Figure  3),  which  gives  the  change  in  average 
frequency  from  one  r  interval  to  the  next  [7].  The  augmentation  incorporated  in  Totvar  (t,T) 
combines  the  sequential  sampling  function  with  others, which  makes  its  sampling  technique  bizarre, 
but  nevertheless  shown  in  Figure  4. 

Although  the  time-domain  presentation  of  the  action  of  Totvar  on  frequency  residuals  seems 
to  give  little  insight,  we  can  use  the  Fourier  transform  of  these  sampling  functions  to  derive  frequency 
responses  that  perhaps  convey  more  meaning.  The  continuous  analog  version  of  Totvar(r,  T)  can 
be  written  in  the  form 


Totvarc  (r,  T) 


1 

2t2T 


du  y  (u)  hy  (u;  t,  r) 


2 


(10) 
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where  hy  («;  t,T)  is  the  sampling  function  for  (<)  (for  the  extended  x*  (•)  record)  in  terms  of 
y{u)  ~  dx{u)/du,  0  <  u  <  T.  Let  Hy{f\t,T)  =  Jq  hy  (u;  t,  t)  We  can  show  that  the 

expected  value  is  given  in  terms  of  Sy  (/),  the  one-sided  spectral  density  of  y  (t),  by 

£;[Totvaic(r,r)]  ==  /“  Wy  (/;  r,  T)  5y  (/)  d/,  (11) 

Jo 

where 

w,  (/;  T,  T)  =  ^  jf  IB,  (/;  f,  t)P  rft,  (12) 

Wy  (/;  T,  T)  can  therefore  be  regarded  as  the  mean  frequency  response  of  Totvar  as  an  action  on 
y  (t).  Figure  5  shows  this  frequency  response  plotted  against  /r  for  r/T  =  0.1, 0.2, . . . ,  0.5.  The 
frequency  response  of  Allan  variance  (dotted  curve)  is  the  limit  of  Totvar  response  as  r/T  —*  0. 
Also  shown  (dashed  curve)  is  the  frequency  response  of  2  •  Totvarc  (T,  T),  regarded  as  an  estimator 
of  Avar  (T/2)  (see  Section  5  below). 

Totvarc  (t,  T)  has  an  approximate  Allan-like  response;  more  importantly,  Totvar  does  not 
have  the  deep  nulls  encountered  with  Avar(r)  near  r  =  T/2;  consequently  it  has  less  variability 
and  hence  better  confidence  as  indicated  by  an  increased  edf.  We  showed  that  the  estimate  {vis-d,- 
vis  E  [Totvarc  (r,  T)])  is  an  unbiased  estimate  of  Avar(r)  for  WHFM  noise  and  slightly  biased  for 
FLFM  and  RWFM  noise.  This  may  be  somewhat  evidenced  by  noting  the  slight  reduction  in  the 
amplitudes  of  the  main  lobes  of  the  frequency  responses  in  Figure  5. 

5  Properties  of  Totvarc  (T,  T) 

Althou^  we  have  defined  Totvar  for  r  all  the  way  up  to  (W®  -  1)  7U,  we  can  realistically  expect  to 
obtain  meaningful  frequency  stability  results  only  r  <  T/2.  Nevertheless,  we  computed  the  mean 
and  edf  of  Totvarc  (t,T)  for  T/2  <  r  <  T.  For  white  and  random  walk  FM,  tvnce  Totvarc  (T,T) 
is  unbiased  for  Avar  (T/2)  and  is  almost  unbiased  (within  a  factor  of  2/  (31n(2))  =  0.9618)  for 
flicker  FM.  This  result  was  motivated  by  comparing  the  frequency  responses  of  2  •  Totvarc  (T,  T) 
and  Avar  (T/2),  as  seen  in  Figure  5.  Unfortunately,  Totvarc  (T,  T)  has  a  smaller  edf  (1.5, 1.126,  and 
1.029  for  WHFM,  FLFM,  and  RWFM)  than  Totvarc  (T/2,  T)  does,  and  hence  the  mean-squared 
error  (bias^  -I-  variance)  of  2  •  Totvarc  (T,  T)  as  an  estimator  of  Avar  (T/2)  is  greater  than  that  of 
Totvarc  (T/2,  T),  so  we  mi^t  still  prefer  Totvarc  (T/2,  T),  or  perhaps  some  linear  combination  of 
the  two.  This  possibility  has  yet  to  be  investigated,  but  other  work  indicates  that  Total  variance 
coefficients  beyond  r  =  T/2  could  justifiably  be  incorporated  in  the  last  T/2  value  usually  reported. 
Summing  all  the  famifiar  “power-of-2”  r- values  in  a  Total  variance  plot  leads  to  exactly  twice  the 
standard  sample  variance,  much  in  the  same  way  that  integrating  an  estimate  of  a  spectrum  also 
yields  the  sample  variance  [18]. 
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Figure  1:  (a)  Extension  of  a  phase  record  by  reflection  at  both  ends;  (b)  circidar  representation  of 
extended  phase  record 


Figure  2:  Expanded  views  of  the  lower  tails  of  empirical  Totvar  distribution  functions  (lower  curves), 
and  the  corresponding  chi-squared  distribution  functions  (upper  curves) 
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Figure  3:  Usual  sampling  ©ven  by  Avar  for  values  of  m  and  n  when  Nx='^0,  rn  <  {Nx  —  l)/2 


m— (N>-1)/1G  m-(N,-1)/8  rTW(IM,-1VA  m— (N^-IVS  m-Nx-2  m-Nx-1 


O  0.5  1  T  .5  2  2.5  3 

f-c 

Figure  5:  Mean  frequency  responses  of  Avar  and  Totvar  as  operations  on  y  (t) 
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Questions  and  Answers 


CLAUDINE  THOMAS  (BIPM):  If  we  want  to  use  the  total  variance,  we  must  first  be  sure  that  there  is  no 
frequency  drift  in  the  data. 

DAVE  HOWE  (NIST):  Oh,  that  is  not  true.  We  no  longer  have  to  remove  the  drift  or  some  coefficient,  at 
any  rate.  We  never  know  the  drift  in  real  data;  in  simulation  we,  of  course,  can  - 

CLAUDINE  THOMAS:  But,  when  you  circle  the  data,  if  you  have  a  frequency  drift,  of  course,  you  will 
have  an  offset. 

DAVE  HOWE:  Well,  no.  As  a  matter  of  fact  -  let  me  put  that  slide  back  on  here,  that  is  a  very  good 
question.  I  will  put  this  up.  This  is  the  original  data,  this  is  20  simulations  of  random  walk  FM.  So  the  data 
duration  extends  with  Index  1  to  1025.  The  procedure  calls  for  extending  that  at  both  ends.  So,  the  vector 
that  you  actually  work  on  is  going  to  look  like  that. 

Now,  in  the  presence  of  drift,  because  it  is  a  mirrored  reflection,  if  there  is  a  trend  in  the  data  at  the  adjoin, 
when  you  are  adding  or  extending  the  data  set,  then  there  will  not  be  any  discontinuity  by  taking  the  second 
difference.  Yes,  you  have  to  think  about  it  a  little  bit,  but  what  we  do  is  we  actually  reverse  the  data  set  and 
we  also  mirror  it.  So,  even  if  there  is  a  racing  trend  through  that  last  point,  or  the  first  point,  then  there’s  no 
discontinuity  in  looking  at  the  variance,  which  would  be  the  second  difference  in  that. 

CLAUDINE  THOMAS:  Thank  you  for  your  answer.  I  also  have  another  question.  Does  the  use  of  the 
total  variance  suppose  that  your  system  is  more  or  less  stationary;  that  at  the  beginning  and  the  end  of  the 
Tau,  of  the  big  T  lines,  the  statistical  properties  of  the  system  should  be  the  same?  Because,  when  you 
circle,  the  statistical  properties  are  completely  different.  I  am  thinking  about  the  total  time  scale  of  course. 

DAVE  HOWE:  Yes.  First  of  all,  by  construction,  we  have  a  situation  where  the  variance  for  the  original 
data  has  to  be  the  same  as  the  variance  for  the  extended  data,  and  that  again  is  by  construction.  We  have  not 
added  any  additional  data;  we  are  using  the  original  data;  and  constructed  it  in  a  way  that  added  more  terms, 
what  we  call  “surrogate  terms,”  which  will  contribute  to  the  answer  that  you  get  in  the  long  term. 
Nevertheless,  getting  back  to  your  remark  about  circularizing  the  data  - 

CLAUDINE  THOMAS:  My  question  was  you  have  a  span  of  data  - 

DAVE  HOWE:  The  treatment  that  we  apply  does  not  have  any  restriction  regarding  stationarity. 
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